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Perturbative QCD analysis of near-to-planar three-jet events. 
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We present the all-order resummed K ou t distribution as a measure of aplanarity in three-jet events. 



1. Introduction 

In recent years much theoretical and experi- 
mental effort has been devoted to the study of the 
so-called jet-shape variables. Such observables 
are useful to describe the spatial distribution of 
final state particles in high energy hadronic pro- 
cesses, like e + e~ annihilation and DIS. For many 
of these (e.g. thrust and heavy-jet mass B, C- 
parameter Q and the jet Broadenings ||) pertur- 
bative (PT) distributions are already available, 
and it is even possible to estimate the power- 
suppressed corrections due to confinement effects 
(for reviews see (§ , |H) ) . 

For two-jet configurations, the kinematical re- 
gion in which the two jets are "pencil-like" turns 
out to be particularly interesting , since here mul- 
tiple soft gluon radiation effects become essential. 
We meet an analogue of this situation for three- 
jet events if we study near-to-planar configura- 
tions. The variable we choose as a measure of 
event aplanarity is K ont , defined as the sum of 
the moduli of the momenta lying out of the event 
plane. 

The aim of this paper is to present the main 
features of -Kout perturbative distribution, while 
all computational details may be found in . 

The calculation is performed at the "state of 
art" level, which means all-order resummation of 
double- (DL) and single-logarithmic (SL) contri- 
butions due to soft and collinear gluon radiation 
effects, two-loop analysis of gluon radiation prob- 



ability and matching the resummed logarithmic 
expressions with the exact O [aq) results. 

In Section 2 basic definitions and notations 
are introduced, together with a discussion on the 
sources of SL corrections. These are examined 
in detail in Section 3 , which contains the com- 
plete expression for K out distribution to SL accu- 
racy. Gluon radiation effects are included in the 
so-called "radiator" , which is described in Section 
4. In Section 5 is shown the final answer for the 
distribution. We finally draw some conclusions 
and present further developments of our research 
(Section 6). 

2. Definition of the observable 

As it is known, a three-jet event in e + e~ an- 
nihilation is generated by a quark-antiquark pair 
which radiates a hard gluon. For kinematical rea- 
sons these three partons lie in a plane. When ad- 
ditional partons are emitted this is no longer true, 
but it is still possible to define an event-plane and 
to measure the radiation out of this plane. 

First we introduce the thrust and thrust-major 
of the event (Q is the e + e~ centre-of-mass en- 
ergy), 
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We choose to fix the z and y axis along the thrust 
and thrust-major axis respectively. Orientation 
of these two axes is determined by imposing that 
the momentum of the most energetic particle has 
a positive z-component, while the momentum of 
the second most energetic particle has a positive 
y component. Hereafter we attribute ph to the 
right hemisphere (left hemisphere ) if pi lz > 
{Phz < 0). Similarly ph is in the up hemisphere 
(down hemisphere) if ph v > (ph y < 0). The T- 
Tm plane is called the event plane, and the sum 
of the moduli of the out-of-plane momenta gives 
the thrust-minor: 



T m Q = J2 



ph. 
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Such definitions, together with energy- momen- 
tum conservation, force the transverse momen- 
tum with respect to the thrust axis ki = [kx^ky) 
to be conserved separately in the right and left 
hemisphere. Similarly the ir-component of kt 
is conserved independently in the up and down 
hemisphere. 

Near-to-planar three-jet events belong to the 
kincmatical region T ~ Tm 3> T m . At the parton 
level such events can be treated as a hard quark- 
antiquark-gluon system accompanied by soft and 
collinear partons. 

At the Born level no radiation is present, the 
three hard partons are in a plane and, as a con- 
sequence, K out trivially vanishes. 

Hereafter we denote by p\ , P2 and p% the energy 
ordered (pxo > P20 > P30) hard parton momenta. 
There are three possible Born configurations, ac- 
cording whether the momentum of the gluon is 
Pi,P2 or P3. They are named by labelling a con- 
figuration with an index 5 which is set equal to 
the gluon momentum index. In order to clarify 
notations, Figure [l] shows the configuration cor- 
responding to S = 3. 

Secondary parton radiation bring particle mo- 
menta out of the event plane, so that K out no 
longer vanish. We therefore study the "inte- 
grated" K out distribution E, defined by 



das 



dT dTu dK 



C{a s ) Ml(T,T M ] 



dK mn 



c 1 

L 2 


1 V 

1 y Cj 
e 2 

A Pi 






* r - ^ ^ 

J z 

63 


c 3 


p 3 


c 4 



Figure 1. The Born configuration with (5 = 3 
for T = 0.75 and T M = 0.48 . The thrust T 
and thrust major Tm are along the z- and y-axis 
respectively. The four regions Cg in the phase 
space are indicated. 



Here M|(T, Tm) denotes the hard matrix element 
for the production of a quark-antiquark-gluon en- 
semble in the configuration 5 for given values of 
T and Tm- The coefficient function C(a s ), which 
has an expansion in powers of a s (Q), takes into 
account all the effects coming from the hard phase 
space region of emitted partons. 

At the DL level, only gluons which are both soft 
and collinear can contribute to E. An all-order 
resummation of such contributions gives rise to 
the Sudakov exponent 



E(X out ) ~exp \-(2C F + C 
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which is the product of three independent factors, 
one for each hard radiating parton. This simple 
structure reflects the fact that a gluon which is 
collinear to one of the three hard partons is not 
influenced by the presence of the remaining two 
hard emitters. Only soft and non-collinear glu- 
ons are able to explore the inter-jet region and 
see the topological structure of the event. These 
interference effects arc typical of three-jet events, 
and give rise to SL corrections to the above dis- 



tribution. The other sources of SL contributions 
are present also in two-jet shape variables. They 
can be obtained by properly treating hard parton 
recoil, by taking into account hard collinear par- 
ton splitting and by determining the argument of 
the running coupling. 

3. SL corrections 

In this section we describe in detail how the 
SL contributions discussed in the previous section 
modify K out distribution. 

Hard collinear parton decay is taken into ac- 
count by simply replacing the soft matrix element 
with the complete Altarelli-Parisi splitting func- 
tions. 

A two loop analysis of soft gluon radiation is 
needed to determine the argument of the running 
coupling. The result is that, following M, in the 
physical scheme ||, for each emitting dipole a s 
depends on the invariant transverse momentum 
of the emitted gluon with respect to the radiating 
dipole. Furthermore, after the k y integration, one 
can safely replace the argument by twice the k x - 
component, which is dipole independent. 

Exponentiation of hard parton recoil effects is 
possible introducing small recoil momenta q a 
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This allows to split parton phase space into a hard 
and a soft part. The hard part is embedded in the 
hard matrix element, while the soft part can be 
exponentiated together with soft gluon emission 
probability via a multiple Fourier-Mellin trans- 
form. The integrated K ont distribution assumes 
now the form: 



S(^out) 



dv 
2iriv 



a{v)= /[d7d/3]/(/?, 7 )e- K( ^ 7) , 
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where 7 and f3 a are the Fourier variables conju- 
gated to qi y and q ax respectively. 

All the effects of gluon radiation are contained 
in the "radiator" 1Z: 



(8) 



This function is given by the two-loop one gluon 
emission probability w(k) (k = (u>, k) is the gluon 
four- momentum), multiplied by a combination of 
the "sources" uiik) and integrated over the gluon 
phase space. The sources ui(k) indicate how par- 
ticle momenta contribute to K ont - They are one 
for each quadrant, since hard parton recoil ef- 
fects are different in each of this four phase space 
regions. Explicit expressions for w(k), for the 
function I(f3,j) and for the sources ui(k) may 
be found in jFj]. 

We stress that, due to the presence of soft inter- 
jet gluons, we expect the radiator to be a function 
of T and Tm ■ 

4. The PT radiator 

The relevant region for the perturbative contri- 
bution to the radiator is v ^> 1. In this limit, for 
a generic configuration 5, equation ($1) becomes 



K 5 = Ci 5) r(fi 2 ,Ql) + CP r(fl 3 ,Q 2 3 ) 
C{ s) f°° dy 
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where r(/j,,Q 2 ) is the DL function 
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It turns out that even to SL accuracy the ra- 
diator is the sum of three "independent" contri- 
butions, each one corresponding to one of the 
three hard emitting parton and proportional to 
its colour charge Ca \ which equals Cp or Ca 
depending whether the parton a is a quark or a 
gluon. 

The three momentum scales Q a are given by 



Ql = Pj r^ 
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They depend on the event geometry and, in par- 
ticular, each one is proportional to the invariant 
transverse momentum of the emitting parton a 



with respect to the other two. This effect is due 
to the presence of soft, but non-collinear gluons, 
whose radiation is sensitive to the topology of the 
event. The factor g a comes from the hard part of 
the splitting functions: its value is 3/2 if parton 
a is a quark and (3q/2N c if parton a is a gluon. 

Hard parton recoil affects the radiator through 
the scales p, (see || for their definition). 

5. K out distribution 

Using the following simplified notations 



-Kout 
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a s {k x ) Q 2 
R'(D) = (2C F + C A )r'(D) 
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we can give the final result for K out distribution: 



£,5 (-Kout) = e 



= e -^(^out) . 
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The SL function !Fs contains the contribution 
of hard parton recoil. At first order in r' it turns 
out to be 



T s = l- In 2 (2C[ 
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It may seem surprising that the contribution from 
the parton lying on the thrust axis gets doubled 
with respect to the other two, but there is a sim- 
ple argument for this, which exploits momentum 
conservation. When a soft gluon is emitted from 
parton 2 (or 3) both the gluon transverse mo- 
mentum and the hard parton recoil momentum 
lie in the up(down)-region, and momentum con- 
servation is possible. On the contrary, when a 
gluon is emitted from parton 1, if its k t lies in 
the up-region, the recoil momentum lies in the 
down region. Because k x is conserved separately 
in the up and down region, there has to be a re- 
coil contribution from the two remaining partons, 
which explains the result in (|l4]). This prediction 




Figure 2. Comparison of fixed order QCD distri- 
bution with numerical simulations based on the 
exact four parton matrix element for 7/3 > 0.1 in 
units of (a s /(27r)) 2 . 



has been found in agreement with numerical sim- 
ulations obtained with the Monte-Carlo program 
EVENT2 0, as shown in Figure 2. 

All effects of soft inter-jet gluon radiation and 
of hard parton collinear splitting are embodied in 
the scales Q Q , which depend on T and Tm- 

6. Conclusions and outlook 

The most interesting feature of K on t distribu- 
tion in ( |l3| ) is its geometry (i.e. T and Tm) de- 
pendence, which allows to measure the effects of 
large angle gluon radiation. 

In order to compare theory with experiments, 
one has to sum over all possible configurations 
and to integrate over T and Tm varying in an 
appropriate range. At this point some problems 
may arise. 

Till now the available data explore all values 
of T and Tm, but the dominant T~l region has 
to be excluded for several reasons. First of all in 
this region our calculation is incomplete because 
the hard matrix element becomes IR singular as 
the emitting gluon becomes soft or collinear to 
the quark-antiquark pair, so that an additional 



resummation should be performed. Furthermore, 
in this region the definition of an event plane be- 
comes problematic, if not even meaningless. Last 
but not least, even if we were able to perform the 
calculation, the result would not be interesting, 
since we lose the richness of three-jet event topol- 
ogy. In order to avoid this problem one should 
look for data in a "restricted" region, e.g. T<0.9. 
From an experimental point of view it is more 
convenient to put a cut in the three-jet resolu- 
tion parameter. If one uses the Durham recom- 
bination algorithm jnj, a suitable choice is, for 
example, y§ > 0.1. 

However, before any comparison with data may 
be attempted, some theoretical work has still to 
be made. 

First of all one has to compare the resummed 
result with numerical NLO computation based on 
the exact two loop matrix clement in order to 
check resummation to O (a^ In 2 ) and to to com- 
pute the coefficient function C(a s ) with the best 
achievable accuracy. 

Then one has to match the resummed result 
with these fixed order calculations to obtain an 
improved prediction in the whole phase space con- 
sidered. 

But indeed the greatest effort shall be made to 
compute hadronization effects, the so-called non- 
perturbative (NP) "power-corrections" , which are 
due to the emission of gluons with k t ~ Aqcd- 
We believe that these correction are of 1/Q type, 
enhanced by a factor \nK out /Q, in analogy with 
what happens for the jet Broadenings J[^|. This 
will be the aim of our research in the forthcom- 
ing period, and we hope that a complete (PT 
and NP) resummed K out distribution may appear 
soon. 
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